
Calhoun 

iniutuiiaiul AKliiv« ou tfit Nilvdl Poi($ra{jua(« School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


1971 

A computer simulation of ASW interactions. 

Griffin, Charles Donald, Jr. 

Monterey, California ; Naval Postgraduate School 
http://hdl.handle.net/10945/15619 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


http://www.nps.edu/ljbrary 


CsMwun is the Naval Postgraduate School's public access distal repository for 
research oiateriels and tnstitutjiooal pubftcatiions created by the NPS community. 
Cathoufii is named for Professor of Mathematcs Guy K. CatHiuo, NPS's first 
appointed — and publi^d — scholar^ author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
MontereVr California USA 93943 



A COMPUTER SIMULATION OF ASW INTERACTIONS 


by 

Charles Donald Griffin 

L __ 














AppAovO-d {jOA puhtic AzlzaSQ.; dLit.'ilbuXlon Luitanutzd. 


1 


f' 




















A Computer Simulation of ASW Interactions 


by 


Charles Donald Griffin, Jr. 
Lieutenant Commander, United States Navy 
B.S., United States Naval Academy, 1962 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN OPERATIONS RESEARCH 


from the 

NAVAL POSTGRADUATE SCHOOL 
March 1971 


^ 7 ^' 
C ' 



"postgp®’JAt;. SCHOOI. 

K0I3T^’'XY, CAx.x 


F, 


9 ': 340 


ABSTRACT 

The model presented in this thesis is a computer, simulation model 
of ASW interactions between a formation of high value group ships, 
protected by some screening ships, and some penetrating submarines. 
The model is designed for use as an aid in improving the ability of a 
proposed screening tactic in the detection of a penetrating submarine. 

A systematic procedure to improve a screen's effectiveness against a 
known submarine threat is demonstrated,, and an example problem is 
worked using this procedure. 
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I. INTRODUCTION 


In the present era of a continuing and growing Sino-Soviet sub¬ 
marine threat, the problem, of effective anti-submarine warfare (ASW) 
to protect a high value group of shipping is of considerable interest to 
the Naval Officer. In recent years a great deal of effort has been 
expended in an attempt to develop an optimal solution to the ASW prob¬ 
lem; but, due primarily to the complexities of the environment and the 
resulting inability to quantitize parameters,, no) such optimal solution 
has been achieved. It is not the purpose of this paper to provide 
another solution technique for the ASW problem, but rather to provide 
a computer simulation model of an ASW action and to suggest one pos¬ 
sible use for the model. 

The model presented in this paper is designed for use as an aid in 
the design of ASW screens by providing measures of effectiveness of a 
screen in the detection of penetrating submarines, the number of attacks 
made by penetrating submarines, and the number of successful attacks 
made by a penetrating submarine. Using the computer simulation 
model, the effectiveness of a simple ASW screen against penetrating 
submarines is maximized within the constraint limitations provided. 
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n. DESCRIPTION 


The model is a computer simulation of the ASW environment, 
utilizing an event store procedure. The model is written in FORTRAN 
IV for the IBM 360/67 of the U. S. Naval Postgraduate School. A copy 
of the program is contained in Appendix I. 

Before beginning the model description, it is desirable to define 
a few terms commonly used in the model. 

A replication is defined as one complete ASW action. The action 
begins with the initial courses, speeds, and locations set by the user 
for the ships and submarine and proceeds until termination of the 
action. 

A run is a collection of one or more replications, each of which 
has the same initial values of initial courses, speeds, and locations 
of ships and submarine. The variation among replications in any run 
are caused by the randomly generated occurances within each replica¬ 
tion. 

A high value group is a group of one or more ships to be protected 
by escort ships from submarine attack. 

A screen is defined as a barrier formed by a group of escort ships, 
distributed about the high value group, whose purpose is to protect the 
high value group from submarine attack. In both actuality and the 
model there are two general screen types, the random patrol area 
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screen and the fixed station screen. In the random patrol area screen, 
screening ships are assigned a circular patrol area with a specified 
radius and a center that is fixed at a relative range and bearing from 
the high value group. The screening ships randomly patrol this area 
at some speed greater than that of the high value group. In the fixed 
station screen, the screening ships are assigned a fixed station relative 
to the high value group and maintain this station throughout the 
replication. 

CPA is defined as the closest point of approach of a submarine to 
a screening ship. 

Detection is defined as occurring when a penetrating submarine 
reaches a point within the detection range of a screening ship. 

A formation is defined as the group of ships that make up both the 
high value group and its screen. The center of a formation is referred 
to as ZZ and may or may not be occupied by a ship. 

The m.odel is designed to simulate an ASW environment. The 
scenario consists of a formation, proceeding along a predetermined 
track and being attacked by submarines attempting to penetrate the 
screen in specific places. The entire ASW action takes place suffi¬ 
ciently far in the first quadrant of an X-Y grid to ensure positive 
locations at all times in the model. The ability of the individual 
screen ship to detect the penetrating submarines is determined by a 
detection range that is uniquely determined for each screen ship in 
each replication. If a random patrol area screen is used, the ability 
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of a screen ship to detect a submarine is also influenced by the location 
of the screen ships within their patrol area;’ i;..e..,, at the extremities of 
the area, the screening ship may be several thousand yards closer or 
farther away from a penetrating submarine than, is the center of the 
ship's patrol area. To use the model, the user must supply certain 
data consisting primarily of the initial locations,, courses and speeds 
of all ships, and some cui*ve of detection range vs., probability of 
detection for the individual ship's sensors. The specific data are 
placed on control cards in a standard format and, a. new set of control 
cards is u£;ed in each run. The flow of simulation.log is contained in 
Figure 1. 

There are five events in the model. These are: 

1. A course change by a screening ship. This event, called 
whenever a screening ship reaches the edge of its patrol area, gen¬ 
erates a new course for the screen ship to continue its random patrol¬ 
ling. 

2. A detection by a screening ship. In this event the location of 
the detecting screen ship, detected submarine, and ZZ are noted at 
the time of detection. 

3. A CPA by a screening ship. At CPA,, the time; CPA range; 
screen ship and submarine identitites; and locations of screen ship, 
submarine, and ZZ are stored in a vector for reproduction at the end 
of the run. 
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Figure 1: Simulation Logic Flow 
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4. An attack by the submarine on a high value group ship. The 
time and location of the submarine and attacked ship for each attack 
are noted in this event. 

5. Termination of the replication. This event locates all screen¬ 
ing ships, submarines, and high value group ships at replication ter¬ 
mination. The mode of termination, either by successful submarine 
penetration or by a sufficient number of screen ship detections on.the 
submarine, is indicated. 

The measures of effectiveness for a specified formation computed 
in the model are: 

1. The number of attacks on high value group ships made by 
penetrating submarines. 

2. The number of times the penetrating submarines are detected 
by screen ships. 

3. The number of successful penetrations made by the submarine. 

These measures may be used to measure the effectiveness of a 

screen. The expression "effectiveness of the screen" refers to the 
probability of the screen ships detecting the penetrating submarines 
and, at the same time, detecting the submarine sufficiently far from 
ZZ to prevent attacks on the high value group ships. Clearly, a screen 
may be designed that will maximize the probability of detection of the 
penetrating submarine; but, with the necessity of minimizing the 
number of attacks made on high value group ships as an added con¬ 
straint, the problem of designing an effective screen is complex. 
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With the measures of effectiveness computed in the model, the user may 
evaluate the relative effect of changes in parameters; e.g^, changing 
the radius of the random patrol area on the effectiveness of the screen. 

It should be noted that the model docs not consider the aspect of 
screen ship attacks on detected submarines. There are two reasons 
for this. The first reason is that the inclusion of screeu ship attacks 
on detected submarines is not considered germaine to the measurement 
of screen effectiveness as previously defined. The second reason is 
that the inclusion of screen ship attacks would greatly complicate the 
model, adding to both time and size requirements. 
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UI. LimTATIONS OF THE MODEL 


The limitations of the model fall into the areas of dimension limita¬ 
tions and ship maneuvering limitations. A more complete description 
of these limitations is contained in the following paragraphs. 

A. DIMENSION LIMITATIONS 

The model is limited to certain maximum numbers of ships and 
replications in each run to preclude using excessive computer storage. 
These limitations are; 


1. 

Screen ships 

10 

2. 

Submarines 

10 

3. 

High value group ships 

10 

4. 

Replications 

20 


B, MANEUVERING LIMITATIONS 

1. Screening ships are constrained at all times to remain within 
their assigned screening stations. In the event of a detection of a 
submarine, the screen ships do not attack but remain in their stations. 

2. Submarine maneuvering is confined to determination of a 
course to pass under ZZ. 

3. High value group ships do not maneuver but maintain course 
and speed throughout the replication. 
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IV. ASSUMPTIONS 


In the model certain assumptions have been made to reduce the 
complexity of the simulation. These assumptions are discussed in the 
following paragraphs with respect to screening ships, penetrating 
submarines, and high value group ships. These assumptions are con¬ 
sistent with the philosophy of model design and the measures of effec¬ 
tiveness derived for the model. 

A. SCREENING SHIP ASSUMPTIONS 

1. Detection ranges are assumed to be Boolien in nature; detection 
occurs with probability one when the penetrating submarine is within 
this range, and detection does not occur otherwise.. The detection 
ranges are generated uniquely for each ship in each replication. 

Detection ranges are generated at the beginning of each replication 
using the following procedure: The model user provides some curve 

of detection probability vs. detection range for each sensor that is 
associated with the screen ships. A number is randomly drawn from 
a uniform distribution on [O, 1]; this number is the detection probability 
for the ship in that replication. The range corresponding to this detec¬ 
tion probability for the sensor is the detection range of the ship for that 
replication. 

2. As indicated above, the model is designed making an assumption 
of a uniform distribution on [O, l] for detection probability. If the user 
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feels that detection probabilities are distributed according to some dis¬ 
tribution other than uniform, the model subroutine generating screen 
ship detection ranges must be rewritten by the user-to reflect this 
different distribution. 

3. Screening ships remain within their assigned stations at all 
times in the model. 

4. Screen ships maintain a constant speed throughout each run. 

5. Upon detection of a penetrating submarine,, the detecting ship 
does not attack the submarine. The detection of the submarine is noted, 
and the replication continues. 

6. Once a screening ship has gained contact on a penetrating 
submarine, the interactions of that screen ship and that submarine 
are terminated; i. e. , as the replication continues,, the screen ship 
does not redetect the submarine or calcvilate a new CPA for the sub¬ 
marine with each new random course change of the screen ship in its 
patrol area. 

B. SUBMARINE ASSUMPTIONS 

1. The penetrating submarine conducts no evasive maneuvering 
in the model. This assumption corresponds to actual submarine 
maneuvering practices. The submarine fire control solution for a 
penetration yields a course to pass under ZZ., Elaborate submarine 
maneuvering in a screen penetration is avoided,, if possible, as any 
maneuvering has a detrimental effect on the submarine fire control 
solution. 
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2. The penetrating submarine conducts no evasive maneuvering 
while penetrating the screen, nor does it attack-screen ships. The 
submarine does conduct attacks on high value group ships when arriving 
within a weapon's release range set b-y- the-Tiaen for each run. These 
attacks allow records of the effects of submarine penetrations to be 
maintained. 

C. HIGH VALUE GROUP ASSUIVIPTIONS 

1. The high value group ships do not maneuver in the model but 
maintain the initial course and speeds assigned by the user for each 
run. 

2. High value group ship number 1 is always located at ZZ. If 
the formation for the high -value group does not have a ship at ZZ, high 
value group ship number 1 is a dummy at ZZ, and the maximum high 
value group size is reduced to nine. 
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V. INPUTS 


Inputs to the model are by control cards,, one set for each run. 

The entries on the control cards are all right adjusted; i. e. , the entry 
values set as far to the right as possible in the column width specified. 
The input cards are ordered sequentially as listed below and follow the 
source deck directly. 


A. CARD NUMBER 1 


The entry values of control card 1 are integer values. 


Columns 


Definition 


1-4 Number of screen ships in run 


5-8 Number of submarines in run 


9 -12 Number of high value group ships in run 

13-16 Number of replications in run 

17-20 Number of detections of submarine by screen 


Range 

I^N^IO 

I^N^IO 

I^N^IO 

1 ^ N ;^20 


ships to terminate replication 


l^N^ll 


21-24 Print screen ship detection range for each 
replication Sr average detection range for 


each ship for the run 

25-28 Print out CPA Sr locations of each ship Sr 


0, No; 1, Yes 


submarine in replication 


0, No; 1, Yes 


29-32 Print out locations of screen ship, sub¬ 

marine, Sr ZZ at time of a penetrating sub¬ 


marine detection by a screen ship 


0, No; 1, Yes 
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Col\amns 


Definition 


Range 


33-36 

37-40 


41-44 

45-48 

49-52 

53-56 

57-60 

61-64 


Minimum submarine penetration speed 
Print out locations of screen ships, sub¬ 
marines, k ZZ at termination of each, 
replication 

Print out event calendar for each replication 
Print out screening ship information for each 
replication 

Print out high value group ship information 
for each replication 

Print out submarine information for each 
replication 

Submarine weapon release range (lOO's of 
yards) 

Print out locations of high value group ship &: 
submarine at weapon release for each 
replication 


O^N 


0, No; 1, Yes 
0, No; 1, Yes 

0, No; 1, Yes 

0, No; 1, Yes 

0, No; 1, Yes 


0, No; 1, Yes 


B. CARDS FOR SCREEN SHIPS 

Cards for screen ships follow sequentially the first control card. 
These entries and those of following cards must be non-negative and 
may be integer or floating point. For each screen ship, there is one 
card with the following format: 
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Columns 


Definition 


1 - 8 

Screen ship number 

9 -16 

Initial course (in radians) 

17-24 

Ship speed 

25-32 

Initial X-position 

33-40 

Initial Y-position 

41-48 

Sensor type 

57-64 

Patrol station radius 

65-72 

Initial X-position of center of random patrol area 

73-80 

Initial Y-position of center of random patrol area 


C. CARDS FOR PENETRATING SUBMARINES 

Cards for penetrating submarines follow directly after the cards 
for the screen ships. For each submarine, there is one card with 
the following format: 

Columns Definition 


1 - 8 

Submarine number 

17-24 

Speed ( - minimum penetration speed) 

25-32 

Initial X-position 

33-40 

Initial Y-position 


D. CARDS FOR HIGH VALUE GROUP SHIPS 

Cards for the high value group ships follow directly after the 
cards for the penetrating submarines. For each high value group 
ship, there is one card with the following format: 
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Columns 


Definition 


1 - 8 

High value group ship numbei* 

9 -16 

Course 


17-24 Speed 


25-32 

Initial X-position 

33-40 

Initial Y-position 


19 








VI. OUTPUTS 


The outputs of the model are of two types,, general statistical out¬ 
puts and specific format outputs. The generaT statistical outputs are 
always provided by the model, while the specific format outputs are 
provided to the user only if they are desired-. Ta have specific format 
outputs, the user codes control card 1 in the appropriate column with 
a 1 for each specific format output desired and a 0 for those specific 
format outputs not desired. EIxamples of general statistical and 
specific format outputs are included in Figures 2,. 3, and 4. 

A. GENET.AL'STATISTICAL OUTPUTS 

1. The number of attacks on high value group ships by sub¬ 
marines for each replication and the total for the run. 

2. The number of detections made on the penetrating submarine 
by scr eening ships during each replication. The upper limit of this 
value is set on control card one for each run. 

3. The number of successful penetrations to ZZ made by penetra¬ 
ting submarines in each replication and the total for the run. 

4. Average detection range of the penetrating submarine from ZZ 
for each replication; average detection range and its standard deviation 
of the penetrating submarine from ZZ for each run. 
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B. SPECIFIC FORMAT OUTPUTS 


1. Individual screen shlp^s generated detection ranges for each 
replication and the average generated detection range for the run. 

2. CPA of each screen ship to a submarine in each replication. 

3. Locations of ZZ, screen ship making detection, and the 
detected submarine at time of detection in each, replication. 

4. Locations of ZZ, screen ships, and submarines at the ter¬ 
mination of each replication. 

5. Event calendar for each replication.. 

6. Screen ship information at termination of each replication. 
This information consists of: 

a. Screen ship number 

b. Course (in radians) 

c. Speed 

d. Ship X-position 

e. Ship Y_position 

f. Sensor type 

g. Generated detection range 

h. Patrol station radius 

i. X-position of center of random patrol.area 

j. Y-position of center of random patrol area 

k. Time of next random course change 

l. Detection status of ship (1 if ship made detection, 

0 otherwise) 

m. Time of last position update 
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7. High value group ship information at termination of each 


replication. This information consists of: 

a. High value group ship number 

b. Course (in radians) 

c. Speed 

d. X_position 

e. Y-position 

f. Time of last position update 

g. Time first submarine passes under ZZ 

8. Submarine information at termination of each replication. 
This information consists of: 

a. Submarine number 

b. Course (in radians) 

c. Speed 

d. X_position 

e. Y_position 

f. Time of last position update 

g. Number of times the submarine has been detected by 
different screen ships. The upper limit of this value is set on control 
card 1 for each run. 

9. Locations of submarine and target ship at weapon release in 
each replication. 
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1. Examples of General Statistical Outputs 

a. Output for each replication 

WEPLAUNCH DETECTIONS PENETRATION 

1 1 L 

b. Output for each run 

PENETRATIONS WEPLAUNCH AVDET STDDEV 

8 12 6057.6 2386.65 

AVERAGE RANGE OF SUBMARINE DETECTION FROM ZZ 
IN EACH REPLICATION 

7765.18 8809.04 0.0 4286.71 8633..7? 6465.39 

2. Examples of Specific Format Outputs 

a. Individual screen ship's generated detection ranges for each 
replication &: the average generated detection range for the run 
SHIP AVRANGE REP 1 REP 2 REP 3 REP 4 REP 5 

1 3049.9 0.0 1097.8 0.0 3380.6 7427.6 

b. CPA of each screen ship to a submarine in a replication 

SHIP SUB TIME CPA SHIPXLOC SHIPYLOC SUB XLOC 
3.0 1.0 73. 7245. 44238. 68293. 37006. 

SUB YLOC ZZ X LOC ZZ Y LOC 
67857. 40000. 64319. 

c. Location of ZZ, screen ship &: detected submarine at time of 
detection 

DETECTION TIME 64 BY DD 4 ON SUBMARINE 1 AT 
X=36733.00 Y=62409.00 

Figure 2: Examples of Outputs Available in the Model 
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d. Locations of ZZ, screen ships, & submarines at replication 
termination 


TIME 

SCREEN SHIP 

X-POSIT 

Y_POSIT 

86 

I 

42180.0 

7Z546.0 

TIME 

MAIN BODY SHIP 

X-POSIT 

Y-POSIT 

86 

I 

40000.0 

68652.0 

TIME 

SUBMARINE 

X-POSIT 

Y_POSIT 

86 

I 

39910.0 

68724.0 

Event 

calendar for each replication 


TIME 

EVENT NO. 

SHIP SUB 

MAIN BODY 

86 

9 

1 1 

1 


f. Screen ship information at termination of each replication 
SmP COURSE SPD X-COORD Y-COORD SONAR 
4 1.62 15. 35434. 72029. 0. 


DETRANGE 

PATRAD 

PATX 

PATY 

CCTIIvlE DS 

4691. 3 

1800. 

35830. 

71074. 

90. 1. 


XYTIME 

86 . 

g. High value group ship information at the termination of each 
replication 

SHIP COURSE SPD X-COORD Y-COORD XYTIME 
1 0.0 10. 40000. 68652. 86.0 

EVT9TIME 

86 . 

Figure 3; Examples of Outputs Available in the Model 
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h. Submarine information at the termination of each replication 


SHIP 

COURSE 

SPD 

X-COORD 

Y-COORD 

XYTIME 

1. 

1.28 

7. 0 

39910. 

68724. 

86. 0 


TIMES DETECTED 

1 . 

i. Locations of submarine & target ship at weapons' release in 
each replication 


TIME 

SUB NO 

SUB X 

SUB Y 

TARGET NO 

TAR X 

74 

1 

37230. 

67927. 

1 

40000. 

TAR Y 

LAUNCH RANGE 




64654. 

4000. 





Figure 4: Examples of Outputs Available in the Model 
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Vn. MODEL USE 


As an example of the use of the ASW simulation model, the prob¬ 
lem of detection for a screen of four ships is: optimized against a 
penetrating submarine threat. The use of the term, optimization in 
the context of this section refers to a search: for'an.improved, though 
not necessarily optimal, solution. The screen is specified to be a 
random patrol area screen with the four screen-ships symmetrically 
distributed about the formation base course in patroL areas whose 
centers are equally spaced. This symmetric distribution of screen 
ships will cover a specific arc in front of the high value group. The 
screen ships must be located a distance from ZZ that is at least 
equal to maximum submarine weapons' launch range. The submarine 
threat is known, with the submarine specified to pass between screen 
stations. The size of the main body is unspecified. (See Figure 5) 

A. PROCEDUT-E 

The procedure used for optimization in this thesis is that of 
response surface methodology (RSM). This procedure focuses on 
the class of simulations for which a single output variable exists, 
and the input and output variables are quantitative and approximately 
continuous. 
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— ^ . — Submarine penetration tracks 

r Radius of patrol area 

+ Center of patrol area 

^ ZZ 


Figure 5: Four Ship Random Patrol Area Screen with 
Specified Submarine Penetration Tracks 
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If there exists a functional relationship between the quantitative 
output variable and the quantitative input variables, this relationship 
may be expressed as: 


y = g (xj, . . . ,xj^) 


where: y denotes the quantitative response variable 

fir 

X. denotes the input variable 
g denotes the functional relationship 
and the value of the response variable y lies' on some k + 1 dimen¬ 
sional response surface. 

Assuming this functional relationship g to be continuous and 
differentiable, we will assume that a Ta-ylor series expansion of g 
can be expressed as: 

w k k 


y = g(X) = + :^ a.X. + ^ ^ aj^XjXj 


+ . . . 

In the above notation, y is a measure of screen effectiveness deter¬ 


mined for a screen by an average of the number of submarine detec¬ 
tions over all replications in the runs; and is the value of the i^^ 
input parameter for the screen. 

The problem using the simulation model is that we do not know 
the functional relationship g, yet we wish to maximize g(X) subject 
to some set constraints; i. e. , maximize y subject to maximum and 

minimum values for X.. The model is used to estimate the coefficients 

1 

of the above Taylor series expansion by sampling. That is, the model 
will be used to obtain values for y for specified values of X^; and the 
technique of least squares can be used to estimate the values of the 
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coefficients of the Taylor series expansion. Having an approxinaation 
for the functional form of g(X), it is possible,, from this functional 
form, to estimate the gradient of g and to use this estimated gradient 
to determine the direction to move for fartheir sampling in an effort 
to maximize g(X). 

If the initial samples used to determine the values of the estimators 
of the coefficients of the Taylor series expajision are taken from a 
rather restricted area, it is reasonable to make the assumption that 
g is linear throughout this area. That is, the functional form of g(X) 
may be written: 

k 

y = gW = ^ a.X. 

and that the gradient of g is: 

(a p . . . , 

The sampling procedure is to observe at least k + 1 values of y at 

k + 1 different values of X and to estimate the values of (a,..,ai, . . . , a, ) 

— U ^ k 

by least squares. A large number of simulation runs with different 
values of the X^'s may be made to determine the values of the a^'s by 
least squares. By considering only the largest and smallest values of 
the Xj's within the restricted area, there are 2 experiments that may 
be performed without redundancy to determine the value of the a^'s. 

The minimum number of experiments that will uniquely determine the 
values of the a^'s in k + 1 dimensional space is k + 1. Due to the high 
cost of computer time, the author did not consider it necessary to use 
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more than k + different values of the input variables to determine 

the value of the a.'s. 

1 

To ease the computational complexity of the solution, the user may- 
code each with +1 for its largest valu€^ and _L for its smallest value. 
B-y coding the input variables as +1, the assumed linear function 
relationship g is transformed to a new linear functional relationship 
of the form: 


y = b + 21 b.x. 
0 111 


where x. = _ ^ + 




A 


+ = maximum X- + minimum X,- 


/\ = maximum X. _ minimum X. 

- 1 1 


and + 2aj ^ + 

= aj 

By using k + 1 experiments, a unique solution for the estimated values 
of b^ through b^^ is obtained by the simultaneous solution of k + 1 
linear equations. 

The estimates of b^ through b^ obtained by sampling are then 

used as estimates of the gradient of the response surface. By using 

the gradient as a guide, one expects to proceed to an improved solution 

k 

V = b + w b.2 
o j 1 

where y denotes the predicted response along the path 
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(wb^, wb2, .... wb^) 

where w is a positive constant to indicate we are moving in the direction 
of increasing values of y. That is, we shall observe a new value of y 
at the point (wb^, wb^, . » . wb^^) and hope that the observed value of y 
will be close to the predicted value of y using y above. It should be 
noted that with increasing w, the accuracy of the predicted response 
may become unreliable unless the unknown functional relationship is, 
in fact, linear; and, thus, this path should be followed with caution. 

As long as observed values for y continue to increase along the path, 
one can continue. If the linear approximation for the unknown func¬ 
tional relationship has seemed satisfactory along the path, the appear¬ 
ance of stationary or decreasing values for y as the path is continued 
to be followed indicates that the linear approximation and gradient is 
calculated. The new path indicated by this gradient is followed until 
either the optimal solution is reached, or a constraint is reached; or 
a new linear approximation and gradient must be calculated. 

B. THE PROBLEM 

- For the problem being considered, the screen is specified to be a 
random patrol area screen with the ship located on bearing 300°, 340°, 
020°, and 060° relative to ZZ respectively. The screen center of 
each patrol area is to be no greater than 10, 000 yards from ZZ and 
no less than 4, 000 yards from ZZ. The submarine is to penetrate 
from 320°, 000°, and 040° relative to ZZ. The number of ships in the 
high value group is unspecified. (See Figure 5) 
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The model provides for each run the number of successful pene¬ 
trations made by the submarine. Therefore, using the number of 
unsuccessful penetrations as a measure of screen effectiveness to 
optimize using RSM, one must first consider what va-riables in the 
model affect this measure of effectiveness. The variables that affect 
the output parameter of effectiveness y are; ship speed,, radius of 
random patrol area, the range of the center of the irandom patrol area 
from ZZ, ship detection range and the location of the ship within its 
random patrol area. The first three variables, ship speed, radius of 
patrol area, and the range of the center of the patrol area from ZZ, 
are variables over which the user has direct control in the model. 

Ship detection range and the location of the ship within its random 
patrol area are random within the model and, therefore, not directly 
controlled. 

For this problem, it was decided to limit the number of computer 
runs to eleven, with ten replications in each run, and to optimize the 
screen effectiveness using no more than those runs. Screen ship 
speed was fixed at 15 knots with the speed of the high value group at 
10 knots. Thus, there are two variables to consider in optimizing the 
screen: radius of the patrol area (X^) and the range of the center of 
the patrol area from ZZ (X^). 

For sampling, the initial area considered was bounded by the 
values of 1, 000 yards and 3, 000 yards for X^ and 8,.000 yards and 
10, 000 yards for X 2 . The coded variables are; 
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X 

X = 1 - 2,000 

1 1,000 

X = ^2 - 9,000 

2 1,000 

The initial sample points and their associated responses were: 


Sample point 


^2 

^2 

y 

1 

1,000. 

10,000 

-1 L 

0.533 

2 

1, 000 

8, 000 

-1 -I 

0. 567 

3 

3, 000 

8, 000 

-1 _L 

0. 667 


The solution of the three simultaneous linear equations yielded 
the following linear approximation; 

y = 0.457 + 0.017x^ _ 0.075x2 

Since the largest observed value of the response was 0..667, it 
was decided to choose w such that the predicted response would be in 
the neighborhood of 0.700. The w used was 41, which corresponded 
to y = 0. 697 . We, therefore, sampled y at the point 
(41 (0.017) , 41 (-0.075)) 

The data for the next observation was as follows: 

Sample point w X. X x x y 

4 41 2,700 5925 w(0.017) w(-0.075) 0.667 0.697 

The value of the observed response for this run was 0.667, a value 
less than predicted but still not less than the greatest previously 
observed value. It was decided to continue along this path to see if 
the solution would improve. Accordingly, a value for w of 49, 
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corresponding to a predicted response of 0.747, was chosen. We 
sampled y at the point 

(49 (0,017) , 49 (-0..075)).. 

The data for this next observation was fallows: 

Sample point w y y 

5 49 2800 5325 49(0..0L7) 49{-0.075) 0.767 0.747 

By examination of the differences between pxedicted and response out¬ 
put for the previous runs, it appeared that llae linear approximation 
was not valid. Nevertheless, as long as the response continued to 
increase, the path indicated by our approximation of the gradient of g 
appeared to indicate the direction to an improved solution. Continuing 
along the path, we next selected w = 58, a value corresponding to 
y = 0. 800, and sampled at the point 

(58 (0. 017) , 58 (-0.017)). 

The data for this next observation was as follows: 

Sample point w X^ X^ x^^ x^ y 

6 ' 58 3000 4650 58(0.017) 58(-0.075) 0.667 0.800 

As the value of the response has-decreased, it was necessary to 
return to the point of maximum response and calculate a new linear 
relationship and gradient. Again, three sample points were used with 
the previous highest response point as one of these. These sample 
points and their associated responses were: 
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Sample point 


X 

2 


X 

2 

y 

7 

1800 

5325 

-I 

1 

0. 733 

8 

2800 

5325 

I 

1 

.0.733 

9 

2800 

4325 

I 

_1 

0 . 800 


The solution of the simultaneous linear equations associated with 
these sample points yielded the following linear approximation: 
y = 0.7665 + Ox^ - 0..0335x ... 

It is of interest to note that the value o£ the estimator correspond¬ 
ing to the value of the radius of the patrol! area is zero. This indicates 
that the value of the radius is at its optimal! value and can, therefore, 
no longer contribute to an improved solution. A value of x^^ =0 cor¬ 
responding to Xj = 2, 300 yards was, therefore, assigned in the final 
runs of the model. 

We decided to choose w = 30 such that the predicted response 

A 

would be y = 0. 800. We sampled y at the point 
(30 (0), 30 (-0.0335)). 

The data for this next observation was as follows: 

A 

Sample point X^ x^ • x^ y y 

10 2300 4820 0 30(-0.0335) 0.800 0.800 

It was realized that the predicted response of sample point 10 was not 
an improvement over the best previous value; but the value of X^ was 
approaching the overall minimum constraint of 4, 000 yards, and a 
more cautiovis step w'as taken. 
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Continuing along the path for the final run, a value of w = 74 
corresponding to y = 0. 850 was chosen. We sampled y at the point 
(74 (0), 74 (-0.0335)). 

The data for this final observation was as fbllawa:. 

A 

Sample point X X x x y y 

11 2300 4085 0 74(-0..0335) 0.833 0.850 

With the completion of the eleven runs of the model, the optimal 
screen for the given constraints was found to have values of Xj, and 
X^ of 2,300 yards and 4,085 yards respectively... The optimal response 
for these va.lues of X and X was y = 0. 833; this means that the op- 

J. c* 

timized screen can be expected to detect a penetrating submarine prior 
to the submarine attacking the high value group 83. 3% of the time. 

With additional runs, the response could possibLy be improved; but, 
with only the variable X^ available to vary and with its value so close 
to its lower limit of 4, 000 yards, the improvement would probably be 
marginal. A summary of the above steps is contained in Figure 6. 

In this sample example, a simple screen that had several con¬ 
straints was considered. By using RSM, the screen variables were 
optimized against the threat of submarine penetration. Using the 
same techniques, the model may be used to optimize more complex 
screens with additional variables. 
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Run No 


Run No. 

^1 

^2 

^1 

^2 

A 

y 

Y 

Comments 

1 

1, 000 

10,000 

-1 

1 


0.533 

Value for path of 
steepest ascent are: 

2 

1, 000 

8, 000 

-1 

-1 

— 

0.567 

(2000 + 17W, 9000-75w) 
iin.Xj^,X 2 space 

3 

3, 000 

8, 000 

1 

-1 

-- 

0. 667 

y=0.457+w(0.005914) 

4 

2, 700 

5,925 

-- 


0.697 

0.667 

w = 41 


2, 800 

5, 325 



0.747 

0.767 

w = 49 point of 
maximum response 

6 

3, 000 

4, 650 


-- 

0.80 

0.667 

w = 58 

7 

1, 800 

5, 325 

.1 

1 


0.733 


8 

2, 800 

5, 325 

1 

1 


0.733 

y=0. 7665+w(0.00112225 
b =0.0 in this solution 

9 

2, 800 

4, 325 

1 

-1 


0. 80 


10 

2, 300 

4, 820 



0.8 

0.8 

w = 30 

11 

2, 300 

4, 085 

__ 

__ 

0, 85 

0. 833 

w = 74 


Figure 6: Summary of the Search for Improved Solution 


37 
































































































































































































VIII. CONCLUSIONS AND RECOMMENDATIONS 


The model presented in this thesis is a simulatton. model of ASW 
interactions between a formation and penetrating submarines. A pro¬ 
cedure has been demonstrated in which the model was used to improve 
the effectiveness of a proposed screening tactic against a known sub¬ 
marine threat. 

An attempt to use the model to measure the effectiveness of 
various screen tactics was examined. The procedure used in this 
attempt was to determine the average detection range of penetrating 
submarines from ZZ for each screen and to compare these average 
ranges by the analysis of variance techniques. This procedure to rank 
screen effectiveness in order of average detection range from ZZ, 
with the screen having the greatest average detection range being the 
most effective, had intuitive appeal. Unfortunately, the standard 
deviations of the average detection ranges were so great, on the order 
of two to three times the differences between average detection ranges 
of the screens, that no statistical interferences could be made with 
this procedure. 

In an effort to use the same general procedure, an attempt was 
made to compare the effectiveness of proposed screen tactics by com¬ 
paring the average time until first detection of a penetrating submarine. 
This procedure also yielded no statistically significant results for 
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the same reason as the previous proposal. The results of these 
attempts are not included in this thesis. 

The model has not been successfully used to directly compare 
the effectiveness of various proposed screening; tactiics-;.but it may 
be used as shown in the preceding chapter to maximize the effective¬ 
ness of a screen, with certain constraints, againsta specific sub¬ 
marine threat. 
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COMPUTER PROGRAM 
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IF(DX1 .GE.0.0.AND.DYl.GE.0.0)GO 
IF(0X1.LT.0.0.AND.OYl.GE.0.0)GO 
BFTA=ATAN(0X1/0Y1 ) 
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33 RPLCSE=PI+THETA 

IF(BETA,GE.0.0.AMD.BETA.LE.(PI/2.0) )Gn TO 80 
IF(BETA.GE.PI.AND.BETA.LE.( {3.0'''^PI )/2.0) )G0 TO 35 
IF(BETA.GT.(PI/2.0).AN0.3 ETA.lt.PI.AND.RELCSC.LT.ANG l)G0 TO 35 
IF(BETA.GT.( (B.O’t'PI )/2.0) . AND. B ET A. L T . 0.0. AND. RELCS E .GT . ANG2 ) GO TO 
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113 RELCSE=PI+TH6TA 

IF(BETA.GE,0.0.AND.BETA.LE.(PI/2.0))G0 TO 115 
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302 IF(MSTAT(7).EQ.O)GO TO 303 
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EA3=ECAL(1,5) 

ECALd , 1 >=999999999 

RETURN 

END 
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LOCATES ZZ (HIGH VALUE GROUP SHIP 1) 





cn 



—. 4 ^ 





< 



CNjrg 




CVCM 








•• 



-- 4 — 

r —4 ^ 




f—H f—H 

0 


fNJCNj 






a: 


•* •* 

alcC 




0:0: 

h- 



<< 




<< 

< 



KK 




»-}- 

a. 


>> 








LOOO 

xoo 




z:</) 

CL 



^(D 





►—1 


xc >0 

OOO 




i /)0 

X 


i-H CD 

ww 




WW' 

1 /) 


000 








4 W 4 4 W 






X 



^ 4 ^ 





UJ 



mm 




mm 

LU 








a: 


mm 






l: 


•• •• 

WW' 




ww 

00 


»—* t—4 

alol 




ala: 




<< 




<< 

0 


>> 

h-h- 




H” h~ 



0000 




4 — 

WW' 

< 


W 4 W 4 







•^0 


^ 4 ^ 



o 


00 

m • 


00 




00 

CL 

^0 

00 

• • 



f-H 

• • 0 

»—1 

^ II 

• • 

CO 

0 

00 


00 0 0 

X 

r—4 

00 

vOsO 

mo 

LU 

cc 

nOvO OOCSj 

00 

wUJ 

sOsO 


(MCNJ 

X 

< 

r—< 


>x 


4 <-^ 4 ^ 

0 f—< 

»—• 


^ 4 ^ 1—4 

X 

L 0 i-< 

4 --.* 

00 


cc 

1 

00 1- oc 

uu 

1 K 

00 h- 

• • 

0 

< 


• • UJOl—H 

LU 


• • LU 

00 

OK 

X 

H 

OC -^ 1 - 

QC 

KO 

00 -- 


K 

CD 

UJ 

00 mini— CO 

0 

UJ • 

00 K 

4 ^ 00 ^—<r-^ 

0 

0 

W 

00 - xxceo 

t /0 

-Jw^O 

0 0 ^ r —4 ^ 4 ^ 

— c\jc\JOO<t'>- 

CO -J 

00 


(vrsjif'.ir.ooOC' - 


• 

C\J C\) r —4 rH X ^ O 

O-^- ^i-OOXO 



< 4 ^ 

— rr OOCC-I — cc 



■!{• -;}■ ocoo—lO'.-t.-t-H + + 00 

4^0 

-J 

0 >tLnLnu^ 4 - 4 * • • 

-J 

-- II O-J 

-- 4-4 4 - + OOU. ► 

^LiJ UJ . .— 4 f—^ 1 

0 * 4 — II 

—J 


•^uJUJLnLn II II II •ujO 

< 

h-i • 

vtiriLlJUJ.-Hr-._J_J II '—■ 

^s. 5 ;<cci—III 


< 

U_.-H.-t STRICT-—.C^LUUJ 


^ U- 

*. •■ 2 : 5 :<cD II II — 

WH-I»-,<C0 II — UJOC'-' 


w w 

WI— iw^vfLT^sCU-l • •ID 

00 

Kr---uj 

H-t»—H-f 4 ^ w fT) ;> > j—•l ““0 

<y0 

II Qc:cit:i— 1 - II II 4 . 

UJ 

<-H II s: 

II II »ttr.r- 4 ooooww II II O'.-tz: 

• c\)<r\) 

UJ 

irx<t^^ir'trt,-.— <— (Srcc*-' 

K 

h— 0 

> >-w-. W 4 r-~* r —1 »» »*W 4 


X<KO 

K 

uj 1— K X ^ ^ w Lu ir uj!— 

< 

(/)^LUK 00(/)w.^X ' II 

II II II X >-►-■•-* l—UJ Or 00-H 

< 

^ II 

II II II oocr:a^D^— 

0 

X 

II II II II 


0 

►-Hir,LnLnLnco<i<:<iiJLU- llo 

C 

II c>-u---'-*-..-tao:>><rcQOocc>>cu_u. 11 0 

0 


|-<CC)OC_J. 


O _J •-'< CCO Cl _J—loot/) t/XOCOQ—i—loo 00 O'-* •-'-JO —I 


0 

00 

K 

m 

0 

r —4 


(\j 

r~i 

0 

0 

0 



r— 4 

1—4 


o 


ooo 


ooo 


1019 TIME2=(FL0AT(ET)-SSV(J,6)) 




I- 













cu 





^4 











X 





a. 



CO 













x" 



• X 





CO 



<f UJ 








UJUJ 





a. 



ococ 





X 



<0 





0 



CO 





cc 



-J- 





0 



ou 








a:*-* 





ill 



KU- 





X 



<ti-< 





-J 



XCD 


CM CM 



<c • 



UJ 




r\jc\j 

>UJ 


CM CM 

00 Q. 


r— 1 1—^ 


•• ► 

2: 


•• »» 

KCO 


<< 



Xi-^ 


_J_J 

»—< 


UJUJ 



OK 



< 




>> 



czci: 

Uu 


>> 


t/)tO 

XO 


<rc 

ox 


coco 


l/)O 0 

a: 


KK 

UJ 





KCQ 



c^x 


XCA 



ClZ) 


Xto 

XI 15 


H -40 


•--0 

tuco 


K-O 

1 — 


coo 


C/OU 

0 


CO CD 

xa 




'W'W 

XUJ 



LUX 


^ -«• 



tux 



OCO 





K 



0 • 





CO 



0 X 


mm 


mm 

cu:z 


mm 

XCOX 




•* 



^ •• 

<C»-^ 4—1 


f—< r—1 



X 


O-J 

K 


<< 



oo: 


W^w 

CLUJC 


UJUJ 


>> 

UJ 


QC'QZ 

*-*zo 




001/) 

Q_i-h 


<< 

Xt-<Q^ 


>> 


00 to 

X-J 


KK 

ool—CO 


coco 



Co: 


^ w 

cc 

m 

'W'w 

o 


a< 



2 OI /0 

1-4 


• 


ouj 

nO • 


ujo:: 



o 

00 


•-C 

CO 

LU CD LU 


oc 

II 

• • 0 

U^'Q 

-J II 

• • 

a: 3 X 

< 

• • 

(N 

00 >to 

XUJ 

WfO 

00 

OOOK 

UJ 

00 

UJ 

00 r-^ 

^K 

C^LU 

nOsO 

l/) 


vOvO 

2: 

(\J 

<< 

<x 


LUO 

> 


•-H 

^ 0 

>0 

K^-h 


OXH- 

CO 


K 

00 K Q ^ 

0 

1 K 

00 K 

•-K 

1 

00 

— 

• • UJ OK 

X^ 


• • UJ 

U. Q 

—* 

• • 

0 

CO 1 - 0 

0 

KO 

00 w 

»-.u_LU 

K 

C^ C^ 

• 

00 (\J(NJK 0 K 


UJ • 

00 mmK 

OOUO 

UJ 

00 x> 

0 

OO—^XX 00 

Xiu 

N 

( 

0 

oc--^xx 

LU «/) 

'W' 

OC-^-^OO 

• 

CM(NJfM(\JOOC C'~ C 

LU 

»)— • 

CMCMmmcDOO 

D.Z< 

K — 

CM CM U J 0 


^-;;- ^ 5 -oc.cc-J ^c c 

-JX 

CMCK 

-X- OCOO-J 

LOOCL 

<X>t^rc OCX-D 

_jNtin(\J(\J + + _J_JU_ O' 

•J 

!l C-J'J-LomrD + + -J—iU- 0 

1—4 

0 ^ 

►f—^ r—1 + 4' 11 11 


♦'Luu.'caca II II II •LurncM 

<to 

-J^ • 

► ►Luujrom II 11 II 0 

<KiO 

_J ^ ^ U .1 LU r-4 1—^— 

<M“ 5 “ 55 ::S<CD-~^-~UJ 02 -' ♦ 

c 

LL m _j _J 2 :5: < ix)-~^^uj >s-UJ 

CD^ 

U.<< 02 :<CD'tlD 

UJ^ 

^ nJ* in sO ZD UJ • K K 

cox 

r— < w' Uj 


COX 


UJH-41—4WW' ^ ^ 

2:> 

^KK 1 ) 11 •m^—J 

UJ 

C\i II SI 

qcoJkK II 11 •“X X 

ujcoo: 

II -- 

'^-'K K 11 11 ^ 


K X 

OroH-(<<w^r'ir',_i_j_jH-<c>-< 

K LU 

^ < 1 : <: 

(—t/”) 00 —X t—'jj lij 

<U 

r-( u_' 1 —r-i—X X 1 -H H- 

<CO CD 

LU 00 X >-LU LU 

-- II 

II II II oo>>>2:>^w:l*^ 

CD 


ii II 11 II ooa:a:a:z ^ 


X }| 

II II II 

u-fvjcvirocvocooooooouLU. n a. 

OX 

C *—LL rn fT', fo r<'. 0 0 <i < <C C 0 0 

oxx 

r*^ i ™4 1—4 L—.’ C—’ 

•-• < CD 0 C _J _J oo oo 00 0 2 : H- 

-JZ^ 

CZK^^ 

<cqoc:-j-jkkkooo 

-JKX 

(—< CD 0 0 —I —1 <401/) 


0 


CM 






r\j CM 


C\i 

CM ^ 





0 0 


0 

0 





r-J 


r -4 

r -4 





oooo ooooo 


57 



1 



1 



I 

i 





in 


LUO- 



LLJ 



X»-t 




2:aj 



HX 




•-x 



00 




HH 



LU 







OCL 




oo 



3 




O^H 



20 




CO 



CjcC 




oo 


fV(V 

»—»o 


CMCM 


00 LU 


m 

H 


90 ^ r«-* 


00 


CM CM 

OLU 


roro 

corv! 

LUOO 


<r<: 

LUO 


<< 

•• 9- 

x< 


LULJJ 

OOO 


LULU 

1—^ f—1 

HCL 



< 






» 

00> 


QCo: 

C3CqC 

LLOO 


0000 

»—» 


<< 

<< 

Ohh 


0000 

XX 


HH 

HK 




HO 



ww 

2a: 


200 

»—1 


200 


OLU 


—0 

• X 



hhC 

►-4 0Q 


mo 

Ql 


000 

c/)0 

H2 



hhO 


<w*w 


OO 



X^ 



^ -K” 

UJ2 



00 LL 




OO 



»—1 




LU 


c^.m 

QlO • 


rOro 

com 

002 



Olulu 


•*> 

9» »» 

»-H>—1 


CM(M 

OCL2 


rOro 

r-^ f—1 

Xo: 


« 

Qi:oo»—• 


<< 

^ w 

H< 


LULL' 

0 H 


LULU 

ococ 

2 



<xo 



<< 

• GO 


» 

LU 0 


qcq: 

HI- 

UJZ 5 


0000 

32 a: 


<< 

w W 

2oo 


0000 

OLU CD 


HH 



0 


<X 3 




a:o 

9* 


>Xoo 

•» 


oc 


CVJ 



ro 


• t 

5:ll 

< 

CO 

X>LU 

< 

OO 

CO 

CO*—* 

OLU 

• • 

oox 

OLU 

• • 

nC O 

Z)0 

0^ 

CO 

^2H 

0^ 

OO 


L/OLU 

ro> 

sQO ^ 

xc 

Nto: 

>OvO -- 

H ---- 

Ql 

r-HOO 

'V.V* }— 

0 

r-I< 

V. H 

UJ oo 

OOO 

00 

LU 

OOH 

H 

LU 

«l-^ • • 

H-• 

0 1 

OO —' 

»— 'LU 

C 1 

OO ^ 

H OO >j- 

LJL< 


• • (\IC\IH 

LLOOO 


• • rOroH 

< co---->t>- 

>—1 

H 

OO xx< 

•-3LU 

H 

OC XX 

O ^COrH-^XU 

02 

OLU 

OO ooc 

0 00 

OLU 

OO ooo 

-J^Oc\J(\jOCuO 

LULL.- 



LUOOOO 


CO—'--oao 

LLC^r-^-J^ UOO-J 

cvx 

^1—>ttr. (M(\!(\!f\J_J_JU- 

CL>-'< 

h->J-IT <\1 (M ro (^1 _J _J UL 

II •- + -J II 

00 : 2 : • 

oc 9 ^ 

*- 5 r -> oc II II i| 

m Cl, 

C< 

•4:' 00 II II ii 

r— ♦ { LU LU 1*^ rH )| 

UJ 

•CcMcvrvjrvj + 4.^^^ 

LU 

• Cf^icomro + + 

ro<< 2 ^^<cD --0 

02 

0^< <UJLur\;(\i>tLP.vO 

<200 

0 - J < <! LU LU ro ro >t ITi sO 

rHUJLJJ»--^-H<CQ(7'f-HLU 


LU U-LU LLI2^ < CO ♦> •• •‘•LU 

KH >—1 

LULL LULU 

22 <0 •. •- •"! 

9 *^-^ ^x-'^ii«'r» »* 

002H 


^^ ^cv CM 0 

OOH 


pT, fOfo: 

•-'>>^^^11 II f-'f-iZ’ 

X'OO 

CM 11 >>i-i—II ii<<<z: 

luoq: 

fTi II ocQiKi- II 11 «<: 

<;<i >—• 

H 0 

CCViOOoO'-'^^v'Mr^vLLUjLLH- 

HOLU 

<rro<< 

ww-COcDLULULU' 

LU II II II II X>-UJUJI— 

<00: 

UJUiOOOO~^«^X 

<q::cd 

LUO H H 

^^X>-^^1 

w .-^ r-l r-l —! O O Z 

OLUCQ 

^2 II 

II II II oo>>>:i: 

OCL '2 

II II 

II II ooQ^arct:; 

><COOCQOr>>C 

OooO 

LL»—CMC\JC\m'MCO000000O 

033 

LU»—,cr',rri(rvfriC)0<<L<i 

oo < CD O O —J-J oo 00 o 

-JZOOO 

*-»HCXOO-J-Jl0l/)00O 

0002 

j—- flO L-^ -U *U ( 


UCCILI 


O 

o 

Csl 


o 

o 

ro 


O 

O 


oouoo 


oouoo 


58 











HH 

rH<^ 

•'SICG 



0 



3 


o ►s; 



0 



00 

CL •‘^ 

f—1 



— 



• 

r H •-< 

wa^ 



m 




'--CC — 

a: t-in 



< 



X 

a:, •‘in 

< j 



LU 



m 

< _j 

f—►-'O 

1—^ 


• 




H^^O 




0 



•» 

1— ♦^XO 

i<:—.cccv) 

Ol 


LU 



• 

X —CCC\) 


< 


• 



C 


OL r—l J 

0 


— 



3 

QL r-+ •‘O 

s: ►_!< 



in 




3 •-0< 

LUO*—O 

LU 


•* 



H 

LUjO*“<CJ 

H-'-titiUJ 

_J 


♦-H 



3 

I—r-HXLU 

•-^CCQ 

< 





LU 

•-^CDO 


u 


0 



> 

OD> 3 •‘ 

<00 00 — 

LU 


< 



LU 

<0000 — 

OOO *-0 

h-o 


0 




OOO ••'O 

UJQ_I(\) 



LU 



•* 

LUO-JOO 

O 

UJ< 


• 



X 



>a: 


0 



in 

— XH 

oinoo< 

luq: 


3 



•• 

omoo< 

< r-HLUh— 

3 


< 


— 

* 

< 1 —* LU ^— 

O *-000 

LUO 


• 


• 

LU . — 

0 •‘COO 

LUO •'21 

30 


CNJ 


0 

3 f\J 

LUO •‘3 


HO 


< 


LU 

►—1 

•-.f-H p*"- »- 

ro^uoK 



LU 


0 

H H 

ro winH 

<i> 

31L 


• 


0 

- XX 

0 •-x 

UJOOO 

»‘-C 


0 



•• •‘O 

LUOOC - 

•^OOOCi^-- 



LU 


LU 

X r-H r-H 

•<oc a: 

— 00 •‘CX. *-00 

00 LU 


• 



in 11 •• 

C\) •'Cv — 


h-5: 


— 


cc 

•> ^ CNJ 


•-Lumooc 

3»—• 



ro 

LU 

V, * 

LULn--0 

CO •-i-H<rHC\.' 

LUh 




Q 

V* —* •- 

• l' <3 r-H 

<rH ^O^w 

> 


HH 

0 

3 

inx 

I-H •-O'^-^ 

UJ<CUJ>K 

uuo 



H 

< 

• — ^sO 

<OLL> 

•"Ul r-H < 

3 


0 


0 

or Xr-H •• 

LUr-l •*< 

oo :z^- 

3 


< 

0 

< 

lu^ 11 ro 


<2>OCc/) 

LUOO 


0 

0 

0 

0 

3>C 0 

UJLUOOrMOO 

3< 


OOLU 




LUOOr- 00 

•• •‘O 

LU 


• 

cr^ 

H 

LUX—X 

— Q •- - 

•‘0 — 0 

Lua: 


00 


3 

OO-^OO 

0 ♦‘C — 

<►— «^,-hO ^ 

0:0 


H3 


LU 

< 0 •• •• 


UJ^iC\ #-00 

03 

ro 

< 

c> 

> 

OCC»— 'CO 

•‘CNJ 

•-OOrHO *-0 

H *—< 

0 

0 • 

0> 

i-H LU 

w-K-< 

— OOf— <0 •‘ 

z:uj •-i-ho*— ' 

00 

0 

CD 

Cr^ 


H30 •^ 

3LU •‘rHO 

UJOO’^'^ 

LL 

(\) 

— < 

cr 

C CD 

^^<ZX 

LUQO^rH 

9-r^<X^O 

00 


3LU 


H 3 

LU<OC0 

w —\ < w 

HO'-'CLa::rvj 

H 

0 

LU • 

cr 

>—* 

>3LU •- 

Ko^c-or 

lu»-h>o<^ 

0 

H 

•0 

• 

C 3 

LU- ClNt 

qc:hh>o< 

^:^oo •‘^H 

CLU 


UL-UJ 

LU 

0 or 

r-H rvj ro — ' >—« 

OOXOO •‘3 

UiCD^ —o< 

LUO 

30 

3 • 

0 

— < 

H3<<< •‘X^ •^ 

— lOD ^*—0 

z:Z)2:oLOH 

00 a: 

LUO 


• 

>^ 3 : .-l(Mr<') UJLL'UJIIjUJ ^XvO — X 

o.r?s om 

O0O0<rH •‘OO 

30 


— ro 

— 

u-i-H- K-z:<<< II II (1 II II 'O'd- ►vj-'J- 

0:loO<i-< 

•‘O —o 


HO^ 

(NJ - 

i-H 

xc *■ luujujuju. 0 C 

•>0 

LUHIsotn •- 

UJO 

d:c> 



<0- II 11 11 11 

11 r-<(\irosj-ir\ o- cQO» in 

oil— 3 :oin 

2:Luoo.-^-->t^ 

3< 

00 • 




^ ^ ^ C\JOZ)(\!OC 

Z LU I/O f—> 

»-* <^D. q: 

* 

00 




10 H H H H H —— CO —— 0 

^ <'-'CL 

Kq[:^q:2:< 

HH 

O-LU 

o< 

0 

UJnC'^ » ^ ► 

••— litTXXXX vO'^» 

1— 

OUJ2<lLUl- 

33 

a • 

<0 

< 

^w)— 

1 -+^. - -ww^wl- ^<^1- 

3LU2:<uj 

OOOKHS: 

OLU 

3 

OUJ 

OrH 

w K—1 _l — I — 1 _! Ci: U- < X LU < c 

CCDOHH 

q:lus:c •^ •^ 

q:zd 

-JLU 

LUw 

roLt 

+1—i-s:3-J _J-J_J 

_JiC<r<<<<i::)—sruni—^1— 

ct:uj30 •‘ 

CDH- 2 : •‘>0 

cco 

- * f —1 

1 w 


>-< —•—cif—<<<<:< ;i oocoOf—i—Qi •-•-■q: Ci 

CDH 3 —X 

3^0 — LU_| 

3UJ 

<LL II U-U- CU- 

11 u_c:cu-iu'oooc!i—uJLUuju-iLLiiJjccc- a:co2r 

330—LU 

OO^-hOOX^^- 

00 00 

o*— 

\p^^, 

H*—‘1 

►-<1—3:LLQ::u.'UJUJuju.xQQ00C^a:3:u.a.3u_0UJ 

001 —»oox 

r-4C\JrO 






p-H 

r-HOsJ 




r-H 

(\) 

1—^ ro 

inrOsC 



o 

o 

C\) 


oeo 

ooo 

(MCSJCNJ 


c 

o 

cv 


oooo 


59 




f 






o 

00 



••2 CD 


z: 

1—1 



0 


*-H 

t 



I 1 > < •• 


oo 

o o 



w(X^ 


< 

X K 



a: ►LO 


UJ 

< 



•* 


oc 

• o 





o 

-- o 





e 

1—4 «—>» 



i^^ocrsj 


1—< 

0^ 






(\J 0^ 



a*f-4 ►o 


LL 

+ cr 



2 ♦*-!< 


o 

“0 O' 



UJO »—0 






K^^UJ 


a: 

-J O' 



p^cnOi 


UJ 

< Q> 





o 

O Q> 



< 0000 ^ 


a: 

UJ 0^ 



000 **0 


O 

• • 



UJO^OJ 



or Of 



Q **hh^ 


-J 

UJ UJ 





< 

• • 



^tnoo< 


»—1 




<f-4UJK- 


K 

f—1 r-4 



0 •'ClOO 



•• ^ 



UJO •*>! 


UJ 

1—4 r—4 



•xi —4 *• 


3 

+ + 

.-^Osl 


rO'-^LTM— 


O 


OO 


— ■*<r> 


UJ 

O'- -- 

0^00 


2UJOOO •• 


00 

J -J 



•- •'OOOo:^ 

OsJ 


OsJ < < 

00 


C<^rsj •'CM •xfM 


2 

ou o 

K-K 


<I<t—xw—X tx 

c 

HH 

Okuj^uj 



UJUJin—ioo 

CsJ 


h- •OsJ • 

OO 


^ •* 1—4 fx-4 (\) 


a: 

oo o 

00 


OsJf-4 •xOxx-^'x-^ 

K 

< 

002:02 



<<OUJ>K 

< 

o 

0 <H< 



UJUJf—4 •*<I<I 



• • 

1—41—4 


x^H- 

U) 

UJ 


•• •• 


«-<2>0000 

»* 

-J 

1—4 1—1 ) 

»—1 r—4 


<UJO0i-4O0 •• 


< 


+ + 


UJ •xQ •* •*— 

o 

o 

T—4 r—4 1—4 Qs ifH 

"^"0 


•x_j 

1—t 


+ + + O' + 



P-4 


H 

O' 

-J-J 


UJN^LH ►CM •• 

o 

2: 

0 

<< 

•x 

►CO 1—4 0 ►O 

rH 

UJ 

0 -J-j«jcr-J 

00 

1—4 



> 

OsJ <<<a'< •— 

UJUJ 


UJOO^'^ ► 

o 

UJ 

oooo'o 

oo 


w •xr-4<^0 

CNJ 


0 oIUJUJCMIJ 

• • 


.-i_j^a cicrosi 


UJ 

J— • • •O' • i^—)CL 

00 CLUJK 

•^-J2 



X 

u) clluhg ‘G •"^ 2 : 

00 C_J 0 

-0<UJ 

2:iiroo ►2H- 

<. 

H 

0 </) O-JOUJGUJ -3_JUJ 

< K • • 

wOK 

UJCC v.^O< 

1- 


• • •□.. • W<^— 

Ol 00 ^^ 

-JUJ 11 

>0X000H- 

00 

00 

— 1 Ql 1 00^^^ •^ _!(_> II 

Z Z'-<'-<ir»<:.Q'~ 

uj{jn<^r^ ►oo 

o 

oc 

O' >Z>2'-<r-lr-I^^Lr,<UJ-~ 


0 II ^ 

••O •x<— X Q 

•• 

UJ 

O' U-l •‘LLi ► •* ) II v*** 

1 —4f—4 k—* r— 4 ^-4 LI -1 •* 

UJH3:OiO ► 


o 

• X.-'X:f-4'D"5'r) «. 

Mt i ti II 


2UJOOr-4w:i^ 

a: 

c: 

G< II < II !i II itr.-'u.;Li!UJoci— 

II ••+GGGCD 

►-4 <1^ CL oc: 

< 

o 

UJ2'^-'5:-)_J_l_J + •>+GGG'-t 

<< 

— GGGGG —1 

H'Q:''^cr x< 

K 


• II II <<<-)< ii:”5-)Z2:z 

II 0 11 -<GGrOiiwwSZ2 2 r 

OUJ2<UJf- 

5: 

K • 

2 00 fx- Q. •-< G G G G CNJ w w w- 1 -X- •-. c 

GCGOGGC 

>—4K-OQ^ 

OC^OHHX 


XUJ 

G C CM UJ G G G LLl (N J O. G G 1—!—1—K m C n'. 0 Q ro C. < < 1—1—1—I—G 

CXU-'XC) ► ► 

-j 

UJ^ 

1- s<<g:gz' 

1— 

200222 H-Q 

CQK2 ►>-; 

-j 

>»—4 

GGCc/iGGGGGGGGGGOCCCGCtOCGGOGGGOOOOGZ 

220—LU-J 

>—1 

UJh- 

*-< G Q 2: C H-i KH G ^ C h-G G G G G G Z: Q G CQ K Q Q G G G 0 G G 

oo—oox»-< 

ro 


1—4 



f~4 r\j ro 


0 

Os.»f- 4 f^ 

0 

OsJ 


CM 

CMCMf -4 

0 

0 

OCO 



OsJ 

ro 

mc^om 


oooo o 


60 












THEM 

IME 

ER 








KCD 








o- 2 : 






-«»r- 


ujLur) 








ujXZ 














cn 


Cl 






-S' 








< 


-JUJX 






•«■ 


>-<}— 00 






o 


:2< 






•t 


0^0 






Os] 


Kuj:^ 








<2:< 






X;. 


Xlu 








1—OUJ 








5 : 





X 

Os] 


ooO*-^ 








CLOOH 





< 



•-H_J 





■>> 

00 


X<00 





0 

a: 


00 ' 





•• 



HX 





(\i 

1 








#-» 


UJUJ 





'v- 

(\j 


UJ>00 



cooo 





a:LULU 



l-HKH 



— 


0 0 ^ 





—Os) 

^CD 


00 



ococ 





(/)►— 


^0 

0000 





Q^<00 




0 

OOO 

0-- 


OUJ 


f\I r—1 

0000 

ou 

a: 

^ Os] 

X— 

LLOrCl 



• • • • 


0-- 

O-W- 

4$* 0 



0 

0000 

0 

• 1 

•-W- 

0 • 

00 < 

m 


• • • • 

h- 


>t< 

• 0 

LU_I 

< 


LULUl—K 


wOO 


>to 

OOOLU 

LU 

^OCDO 

COO,U_I 

0 



r-« 

cr^oro 


CTr-lCDO 

• • • • 

0 



Os] +0sJ 


(\J 


r —1 r—4 


OslX 

OsJOO 


o<< • 

c 

r-lr-l •<- —* 

X XX X 

rvj 

-W-w 

^ 0 

Q ^^CMCyO 

ocLxa: 

LU 

<<—*CC I-I 

xxxx — 


+ 

^•«-o 

oc 

o<t 

•> 

LULUCXi ^‘OCG 

xcxxo 

• 

0 

--r-IQC 

XOsJif^ 

::5X 0 

r—< 


• • • • • 

K 

0 <<rM 

<<< 

0 0 ^ 

ujoLu:^ 

< 

>>*«• < X 

cox C c\J 

*0 

\S\ v 

2—- 

0 — 

^OOOLU 

LU 

— 00 00!-<IJJO'X 

22^ 2*^ 

• 


2-^ 1 

^ --000--C 

2a:-J 

»<• 

-- 1 1 <r--Or-l 

<<<<--^ 

»—lO 

C <1—<LU 

<K^LU 

-- OOCK • 

O0<X< 

r —1 

Cvj^^O —< 

• • • • 


K OCX^OO 

OOC^OO 

- •O^QCO 

^a:oo 

•> 

>tLn<iooc< 

CCOCCLO 

20 


--OCLO 

CL^r-l-X" QCOnO 

1 - 0 

K 

• •< 

• • • • ^ • 

CC2 

0 OOOOw^OO 

0000 7$- 00 

+ ^ 1 OJOOO^i* 

<ia: H 

LU 


OOOOOCM 

X< 

0^0 -f OOLU 

--0 + 0LU 

aCL^OoO + 2 



*00 

^ • • • • 

•?r • 

^OOr-iOO 

>— 1 ;^ #0 

►-H 1 X—i^ox 

luluxlu 


CO 1 LULUo: •ox 

C LUf— 1— LUrOi— ' 

r-H 

•-• 1 <0 II 

CLir <OsJ II 

xo •-a: 1 ^Oi 

2XLJ-*> 

►“•rn 

w^^C?OLU< 


K 

CLi— 

1 


LUH2LLJ 

00 ► 

.--)C\l>>00O 

(\1 • • • 9^^ 

CL • 

• a-};- •«• <(\joo-^ -K'—(\io 

• 1 <ooosjir,^ 

0 

C) ^ 00 OO 11 —J r-H 2 C\l X. C\J ‘ r—1 "f* "f* 

•«• 1—r^u-—'Cir\ II •'iTi 

^OLTi 11 *'ir\ujQ-UJ-|}- •>;• • II 

2LULU 


LU w wO cr CO < 

^ < <i<: <>— iH-i 

^C_J 

0 II • •LUr-l 

II • •LUr-l 

O'—^O-K- 0>4- 

H»-'XX 

a. 

^cc II II oucot—C Q-C <<<: <xx X • *0 

• <fMOt/:<c<(NjOoo<rc 

• X> •LU—'LU 

XXh-H 

_j<i >w-i-(.-i_jc/^qd<;H' 


ii cvjoi-osi'-^oLuKs^^ouji—II 5 : 

UJf— 

JK v'/'w <!££:«:'— w ti 

II WW.W.W- II II 


w-s: 11 II 

5: II II L> 0 w. 

W II 

>‘--'LL2 

C II 

II II <ICU0LUU.i-.0 

•-'CU_U-U_U.--'CM 

II LLOLU<f-^001U>O<I »-^OOX>OLLX II II II (XK 

tu:20»-' 

0<CD0< OjoO^-hj-hXOXO^-^^^^I—KXC 3»—'C:*—'C;<a:lQonOO< ccoooc 

•-hCDLUOXO^ 



CvJ 

CSJ CO 


0 

0 

r-l Osi 



0 



CvJ 

0 

0 Os] 



r—1 




m 

Ln 


OOOOO 


61 













>- 









000 


















zo 

0*- 




- 





•-> 






—* 



h-UJ 









Dec 









LJJQ. 









K 






CD 



LUK 






•K* 



00 






< 



Z 






* 



CL 






0 



UJOO 






• 



x< 






CVJ 



h-X 


















LJJQ. 

m 





— ^ — 




< 








LLX 

UJ 





CO f\i 



OOO 






* ^ 




CNJ 





< 



LJJ 2 : 

< 





^ 00 



XUJ 

UJ 





0 or 



»-hLJJ 






t 



\-cc 

f—^ 





CM 1 



0 

< 








LDOO 

UJ 





V, CM 



1 — 

•‘O 





-- * 



<UJ 

H-O 








cx:x • 

UJ 





^CM —CO 



UJh-UJ 









2: 2 


X 





— 


UJU-h-h 


0 


0 


000 0^ 



Ohhcx 

03 

2 


ro 


41 - X C\J 

X — 


<i 

20 

C 

m 



O'- 0^ 

^ 0 




X 

< 

0 


• 1 

0 • 


h--JCD 

a.— 

0 

X 

1 - 


>t< 

• 0 


2 X 






^(\j W'W 

CM >tO 


OOoo 


X 

CM 

C 


+ -«• +'- 

0 '-'O 


LJJ 

-.i •- 

00 

< 

0 


——— 

vO + CM 


^O-ULl 


cc 

X 



CM CO C^JOO 

— *{{' 


4 X 

O'— 

3 


<\i 


C 

0 —— CMOO 


<XK 

--or 

c: 

r*H 

C K 


—+ —-^-o 

K —X 


000 

o< 

0 

< 

• t 

0 


cnosj-Tr ^ 


2 

• K 


X 

rvj K 

c 

<<I(\j <<IC 

a cD-jj-o^ 


002c 

f-H 

1 - 

»• 

V, 0 

nO 


0 '-— 

ITi 

H-,UJ 

• X 

X 

f-H 

• 


2 —2-- 1 

^ —LOO '—0 ID 

HUJ 

UJ 2 

1-^ 

X 


♦-M #-H 0 

c 

<}—<Xi <K^X 

— OCO(— •!—X 

ujs: 

ooro 

-JO. 

2 

m 

— »-H • 

1- 

OX^oo OX»—'00 

*0_ lXOx2 


20’-" 

•»—I 


X 

-- X 20 


--O—'O W. 0 X 0 

XOr-XGT'D + ' 


♦“•OOh- 

—^ • 

X 

2 

' 1 0 o: 2 

c 

OOOO '—00 OOoo-JJ-LO 

1 1 1 OOOO-^^-N^-h- 

LJJ^ 

u. 0 

CNJI— 

c 

H-• 

CL — • 3 < 

0 

--O + 00 X — o + cx 

C^*>+ 2XX 

2 : II 

LUUJ 

1-hO 


1— 

+ • 


OO'-HOO »— 'O*— • *0 

i-hh-«0C_|w.oX2 II 

1 — 

UJXK 

► • 

X 

X 

Cl ^,-h 

»—t 

1 <0 H CL-«- <C\J 11 

XX ^‘X 1 ^ 


21 -UJ 

*—H LO 

2 


»—• CL — H 

X 

^<'— X—r-r 1 

-- 1 -HO—M* '-K-H 


Q 

<UJ 

>—• 

LuoxeviO 1 1 a • 

»a* ■}><■ <CNOG-Ji- •«■ ^'<\;o 

• XCOOCMlD'-'Xi-J 

11 

1—(X 

1 X 15 : 

1- 

20 

.aa*!—ooLL 

'—C UN II •'sO^ClT. II r-vCl—■>;• ♦ II II •“ 

10 *—^ 

XO< 

wK-1 


00 

1 CNJ*-''—c_i 

0 

II • •UJ*— 4 It • •Xf-H 

0 X >—C ■}{• 0 'C^ UN —< 

LJJ< 

Oll 

>K 

X 

coO'.^xx » *0 

• < CNJ C (/^ < C < rvj 0 c/*'< 0 

♦—> .U- — UJU< 

SILL) 

a: uj 

00 «J 

cc 

Xl— 

^ 1 — II ^-'^^cNjC’i—G 

2 ;— >— ou-h-s:— >— oLut—G*—i/^cNj— II s:s:uj 

t w* 

03<O 


0 

X 

II X 11 11 


s; II II oO'— 5 : II II 00 — 


K> 

X^clC 

LLX 

h- 

II xox<r->oox>>o<—^oox>oxm II ii iixi— i—> 

^00 

0002 


00 

OOCD0XI—l-^O>—'O 


G < a; 0 </) 0 0 < a: Q 0 •— CD u. 0 X c _i _J 00 


i-H 

CO 

0 

0 

f—4- 

fO 

CM 


ro 

0 

0 

0 




vO 

vO 

0 


ooooo ooo oooo 


62 


SUBROUTINE FIND IS CALLED TO GENERATE TIME OF EITHER DETECTION OR 
CPA FOR THE SCREEN SHI^> ONLY IF THE SCREEN SHIP HAS NOT PREVIOUSLY 
MADE A DETECTION ON THE SUBMARINE. 
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STORE TIME OF NEXT COURSE CHANGE 
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STORE TIME OF NEXT COURSE CHANGE 
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APPENDIX B: LIST OF PROGRAM VARIABLES 


This appendix contains the name and definitiorL of th-e symbols 
appearing in the computer program. This list contama only those 
variables using storage space in the program and does not contain 
those used only for temporary storage. 
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Name 


Definition 


CPA(I, J, K) 

CPA information for the ship onr. the submarine 

in a replication. The values hr this vector, are: 

CPA(I, J, 1) 

Screen ship number 

CPA(I, J, 2) 

Penetrating submarine number 

CPA(I, J, 3) 

Time of CPA 

CPA(I, J, 4) 

CPA range 

CPA(I, J, 5) 

Screen ship X location at CPA 

CPA(I, J, 6) 

Screen ship Y location at CPA 

CPA(I, J, 7) 

Penetrating submarine X location.at. CPA 

CPA(I, J, 8) 

Penetrating submarine Y location at CPA 

CPA(I, J, 9) 

X location of ZZ at CPA 

CPA(I, J, 10) 

Y location of ZZ at CPA 

DECAL (I, J) 

Dummy event calendar, a copy of all events entered in 

the event calendar. The entry in DECAL contains: 

DECAL(I, 1) 

Event time 

DECAL(I, 2) 

Event number 

DECAL(I, 3) 

Number of screen ship in event 

/ 

DECAL(I, 4) 

Number of penetrating submarine in event 

DECAL(I, 5) 

Number of main body ship in event 

DESKIL 

An index of successful attacks by a penetrating sub¬ 
marine on a main body ship in a replication 
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Name 

Definition 

DSSV(I, J) 

Dummy submarine vector containing, initiailinformation 

DSSVd, 1) 

tVi 

Number of penetrating submarine 

DSSV(I, 3) 

Speed of submarine 

DSSV(I, 4) 

Initial X location of submarine 

DSSV(I, 5) 

Initial Y location of submarine 

DSTAT(I, J, K) 

til 

Detection range from ZZ in the 1 ireplication for the 

R^^ submarine by the screen ship 

DSV(I, J) 

Dummy screen ship vector containing initial'information 

DSV(I, 1) 

th 

Number of D screen ship 

DSV(I, 2) 

Initial screen ship course (in radians) 

DSV(I, 3) 

Screen ship speed 

DSV(I, 4) 

Initial screen ship X location 

DSV(I, 5) 

Initial screen ship Y location 

DSV(I, 6) 

Screen ship sonar type 

DSV(I, 8) 

Screen ship patrol area radius 

DSV(I, 9) 

Initial X location of center of assigned patrol area 

DSV(I, 10) 

Initial Y location of center of assigned patrol area 

DTARd, J) 

Dummy main body ship vector containing initial in¬ 
formation 

DTARd, 1) 

f Vi 

Number of main body ship 

DTAR(I, 2) 

Initial main body ship course (in radians) 

DTAR(I, 3) 

Main body ship speed 

DTAR(I, 4) 

Initial X location of main body ship 

DTARd, 5) 

Initial Y location of main body ship 
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Name 

Definition 

EAl 

Number of a screen sLip 

EA2 

Number of a penetrating submarine 

EA3 

Nximber of a main body sbip 

ECAL(I, J) 

entry in event calendar 

ECAL(I, 1) 

Event time 

ECAL(I, 2) 

Event number 

ECAL(I, 3) 

Number of screen ship in event 

ECAL(I, 4) 

Number of penetrating submarine in. event 

ECAL(I, 5) 

Number of main body ship in event 

EN 

Event number 

ET 

Event time 

KT 

Index of total number of events entered’in DECALi(I, J) 

MAXEV 

Maximum number of events allowed in ECA1 j(I, J) 

MBKIL 

Number of successful submarine penetrations in 

replication 

MBKILL 

Number of successful submarine penetrations in run 

MSTAT(I) 

Control vector for each run 

MSTAT(l) 

Number of screen ships in run 

MSTAT(2) 

Number of submarines in run 

lvlSTAT(3) 

Number of main body ships in run 

MSTAT(4) 

Number of replications in run 

MSTAT(5) 

Number of detections of submarine by screen ships to 

terminate run 


76 







Name 

Definition 

MSTAT(6) 

Print screen ship detection: range for:-each replication 

it average detection range fojr each.s.creen ship during 

run 

MSTAT(7) 

Print out CPA & locations of each ship &: submarine in 

replication 

MSTAT(8) 

Print out locations of screen ship; submarine, and 

ZZ at time of detection of a penetrating submarine by 

a screen ship 

MSTAT(9) 

Minimum submarine penetration speed 

MSTAT(IO) 

Print out locations of screen ships,, submarines, & 

ZZ at termination of each replication 

MSTAT(11.) 

Print out event calendar for each.replication 

MSTAT(12) 

Print out ship vectors for each replication 

MSTAT(13) 

Print out main body vectors for each replication 

MSTAT(14) 

Print out submarine vectors for each replication 

MSTAT(15) 

Submarine weapon release range (lOO's of yards) 

MSTAT(16) 

Print out locations of main body ship Si submarine at 

weapon release for each replication 

MTARK 

Index of submarine attacks on main body ships 

PI 

3.14159 

R 

Index used in subroutine FIND to isolate sections of 

the subroutine after the first call’on the subroutine 
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Name 

Definition 

SONAR(I, J) 

The generated detection range of:the screen ship 

for the replication in a. run. 

SONAV(I) 

The average o£ generated detection.ranges of the 

SSV(I, J) 

screen ship in a run 

th 

Information vector of the B snbmarine in a replication 

SSV(I, 1) 

Number of the submarine 

SSV(I, 2) 

Submarine course (in radians) 

SSV(I, 3) 

Submarine speed 

SSV(I, 4) 

Submarine X location 

SSV(I, 5) 

Submarine Y location 

SSV(I, 6) 

Time X-Y location fixed 

SSV(I, 7) 

Number of times detected by screen ships 

SUBKIL 

Number of times the penetrating submarines detected 

SV(I, J) 

th 

Information vector of the I^"^^ screen ship in a repli¬ 

SV(I, 1) 

cation 

til 

Number of the I screen ship 

SV(I,2) 

Screen ship course (in radians) 

SV(I,3) 

Screen ship speed 

SV(I,4) 

Screen ship X location 

SV(I,5) 

Screen ship Y location 

SV(I,6) 

Screen ship sonar type 

SV(I.7) 

Screen ship generated detection range 

SV(I, 8) 

Patrol area radius 
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Name 

Definition 

SV(I, 9) 

X location of center of assigned patmL area 

SV(I, 10) 

Y location of center of assigned patrol area. 

SV(1.11) 

Time of next course change 

SV(I, 12) 

Number of submarines detected in replication 

SV(I, 13) 

Time X-Y location fixed 

TAR(I, J) 

til 

Information vector for the I main body ship 

TAR(I, 1) 

Number of the main body ship 

TAR(I, 2) 

Main body ship course (in radians) 

TAR(I, 3) 

Main body ship speed 

TAR(I, 4) 

Main body ship X location 

TAR(I, 5) 

Main body ship Y location 

TEMP(I) 

Used to aid in ordering event calendar as a temporary 

storage location for an event 
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